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Abstract. A sample of 32 galaxies with ac curate distance 
moduli from the Cepheid PL-relation (Lanoix 1999 ) has been 
used to study the dynamical behaviour of the Local (Virgo) su- 
percluster. We used analytical Tolman-Bondi (TB) solutions 
for a spherically symmetric density excess embedded in the 
Einstein-deSitter universe (qo ~ 0.5). Using 12 galaxies within 
6 — 30° from the centre we found a mass estimate of 
1.62A/viiiai for the Vir go clu ster. This agrees with the find- 
ing of Teerikorpi et al. (1992) that TB-estimate may be larger 
than virial mass estimate from Tully & Shaya (1984). Our con- 
clusions do not critically depend on our primary choice of 
the global Hq = 57 kms^^ Mpc^^ established from SNe la 
(Lanoix 1999). The remaining galaxies outside Virgo region 
do not disagree with this value. Finally, we also found a TB- 
solution with the Hq and go cited yielding exactly one virial 
mass for the Virgo cluster 



Key words: Cosmology: theory - distance scale - Galaxies: 
distances and redshifts - Galaxies: general - Galaxies: kine- 
matics and dynamics - Local Group 



1. Introduction 

When studying the dynamical properties of a structure like the 
Local (Virgo) supercluster (LSC) one obviously cannot use the 
standard solutions to the Friedman-Robertson- Walker (FRW) 
world model which are valid in a homogeneous environment 
only. General Relativity does, however, provide means for this 
kind of a problem, namely the so-called Tolman-Bondi (TB) 
model. While no longer requiring strict homogeneity, spheri- 
cal symmetry is a necessary requisite. In the present paper we 
address the question how well the velocity field in the LSC fol- 
lows the TB -prediction, i.e. how well it meets the requirement 
of spherical symmetry. 



1992 



Previous studies (Tully & Shaya 1984 ; Teerikorpi et al. 
hereafter Paper I) provided evidence for the expected 



velocity field if the mass of the Virgo cluster itself is roughly 
Send offprint requests to: T. Ekholm 



equal to its virial mass. However, the Tully-Fisher distances 
used caused significant scatter Now we possess a high-quality 
sample of galaxies with accurate distances from Cepheids. 

Tolman ( 1934 ) found the general solution to the Einstein's 
field equations for a spherically symmetric pressure-free dust 
universe in terms of the comoving coordinates. The metric can 
be expressed as: 

i?'(r,r)2 



ds{r, t) 



dr' 



1 



-dr^ ~R{r,TYdn'- 



(1) 



where dVi^ — dO^+sir? 6d(jP and /(r) is some unknown func 
tion of the comoving radius r. R corresponds to the usual con 
cept of distance. Integration of the equation of motion yields: 

F{r) 



R' 



R 



fir). 



(2) 



R refers to dR/dr and R' to dR/dr. F{r) is another arbitrary 
function of ?•. 

Bondi ( 1947 ) interpreted Eq. g as the total energy equation: 
/(r) is proportional to the total energy of the hypersurface of 
radius r and F{r) is proportional to the mass inside r. In this 
interpretation it is clear that one must require F{r) > 0. Eq. g 
integrates into three distinctive solutions. Two of them are ex- 
pressed in terms of a parameter rj, the development angle of the 
mass shell in consideration: 
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for negative energy function f{r) < 0. The third solution is for 

fir) = 0: 
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It is important to note that the solution has three undefined 
functions: f{r), F{r) and To(r). One of these can be fixed by 
defining some arbitrary transformation of the comoving coor- 
dinate r. In our application we set tq (r) = and interpret r as 
time elapsed since the Big Bang, i.e. we equal t with the age 
of the Universe Tq. 

To find a TB-prediction for the observed velocity, we use 
the formulat ion of Paper I, based on the solutions given by Ol- 
son & Silk (1979). In this approach the development angle ry 
needed for the velocity is solved with the help of a function: 



A{R,To)^ 



^GM{R)To 

i?3/2 



(6) 



where R is the distance of spherical mass shell from the ori- 
gin of the TB-metric, M{R) is the mass contained within the 
shell and G is the gravitational constant, ry is solved either 
from Eq. 7 or Eq. 8 of Paper I depending on value A{R, Tq). 
Ekholm (1996; hereafter E96) calculated in his Appendix B 
the exact value of A dividing the family of solutions into open, 
hyperbolic solutions (/(r) > 0) and closed, elliptic solutions 
(fir) < 0): A = ^/2/3 « 0.47. 

2. Details of the model used 



In Paper I, Ekholm & Teerikorpi (1994 hereafter ET94) and 
E96 one chose as a starting point a "known" mass M{R) for 
each R inferred from some density distribution given as a sum 
of an excess density and the uniform cosmological background. 
We assume: 



p{R) 



3gggo 
AttG 



{l + kR-"). 



(7) 



Hq is the Hubble constant, go is the deceleration parameter 
k and a define the details of the density model and will be 
explained shortly. 

E96 further developed this formalism by defining a "local 
Hubble constant" Hq = Kosmlcf =1), where d is calculated 
from Eq. |l^ By setting d = 1 we can fix the mass excess 
k' independently of the density gradient a. In this manner the 
quantity A can be expressed in terms of d and qq (Eq. 7 in E96): 



A{d,qo) - Ciqo)y/qo{l + k'd-'-). 



(8) 



The factor C{qo) depends on the cosmological FRW world 
model chosen (cf. Eqs. 8-10 in E96). The mass contained by 
a shell of radius d is (Eq. 4 in E96): 



M(d) = :^goi?Virgorf'(l + k'd^'^). 



k' (k normalized to the Virgo distance) provides the mass ex- 
cess within a sphere of radius d — I and can be fixed from our 
knowledge of the observed velocity of the centre of the mass 
structure and of the infall velocity of the Local Group (LG) as 
follows. Let V^Jjgo be the observed velocity of Virgo and V^q 
the infall velocity of the LG. The predicted velocity of a galaxy 



with respect to the centre of the structure in consideration was 
given by E96 (his Eq. 1 1): 
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7 = (K?So + ^liG)/VV°h. 
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(11) 



(?!)(?7o) is the angular part of Eqs. 15 and 16 given by ET94, 
where r/o is the development angle corresponding to the given 
value of A. Hence for each V^q one may fix k' by requiring: 



"(1) = vvTgo 



(12) 



One clear advantage of this formulation is - being perhaps oth- 
erwise idealistic - that the value of k' depends only on the two 
velocities given above for a fixed FRW world model {Hq, qo). 
Thus a will alter only the distribution of dynamical matter (the 
larger the a the more concentrated the structure). 

Finally, we need the predicted counterpart of the observed 
velocity of a galaxy. In the present application it is solved from: 



T^prcd(dgal) - ^V^rgo COS 9 ± v{d)^ 1 - sili^ Q / d^ , (13) 

where dgai = -Rgai/^virgo is the relative distance of a galaxy 
from the LG, O is the corresponding angular distance and d, the 
distance to the galaxy measured from the centre of the struc- 
ture, is evaluated from: 



1 



C 



2d<ral COS 6. 



(14) 



The (— ) -sign is valid for points closer than the tangential point 
rfgai < cos Q and (+) -sign for dgai > cos 9. 

3. The sample 

Eq. |l^ reveals one significant difficulty in an analysis of this 
kind. In order to find the relative distances dgai one needs good 
estimates for the distances to the galaxies i?gai. Such are dif- 
ficult to obtain. Distances inferred from photometric data via, 
say, the TuUy -Fisher relation are hampered by large scatter and 
as a result, by the Malmquist bias which is difficult to correct 
for Distances inferred from velocities using the Hubble law are 
obviously quite unsuitable. 

The HIPPARCOS satellite has provided a sample of galac- 



tic Cepheids from which Lanoix et al. (1999a) obtained a new 
calibration of the PL-relation in both V and I band. Lanoix 



( 1999) extracted a sample of 32 galaxies from the Extragalactic 



(9) Cepheid Database (Lanoix et al. 1999b) and, by taking into ac 



count the incompleteness bias in the extragalactic PL-relation 



(Lanoix et al. 1999c), he inferred the distance moduli for these 



galaxies with 23 based on HST measurements and the rest on 
groundbased measurements. It is also important to note that the 
sample is very homogeneous and the method for calculating the 
distance modulus is quite accurate. The photometric distances 
found in this way are far more accurate and of better quality 
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than those derived using Tully-Fisher relation, a fact compen- 
sating the smallness of our sample. 

We intend to use these galaxies together with the TB- 
model described in the previous two sections to study the gross 
features of the dynamical structure of the LSC. Because we 
need only the observed velocities, T4bs and the Cepheid dis- 
tance moduli /i, we can avoid the usual difficulties and caveats 
of other photometric distance determinations. The velocities 
were extracted from the Lyon-Meudon Extragalactic Database 
LEDA. By observed velocity we refer to the mean heliocentric 
velocity corrected to the LG centroid according to Yahil et al. 
(|1977|). 



Name 



e 



Vol 



d2i 



Vprcd, 1 d20.7 Vprcd, 2 



We need, however, some additional information: V^JJ^q, 
Vj^Q and i?viigo- For the observed velocity of the centre of the 
LSC the value preferred by Paper I is V^f^^^ = 980.0 kms^^. 
For the infall velocity of the LG into the centre we choose 
= 220 kms~i (Tammann & Sandage |1985|). It is worth 



K 



LG 



noting that though this value has fluctuated, this relatively old 
value is still qu ite co mpatible with recent re-evaluations (cf. 
Federspiel et al. 1998[ ). Furthermore, the fluctuations have not 
been very significant. We also need an estimate for the distance 
of the centre of the LSC, i?virgo- One possibility is to estab- 
lish the Hubb le constant Hq from some independent method. 
Lanoix ( |1999| ) found using SNe Ia's:[| 



Ho ==57±3kms~^Mpc^ 



(15) 



This global value of Hq is in good agreement with our more 
local results using both the d irect and inver se Tul ly-Fisher re- 
lations (Theureau et al. |1997|; Ekholm et al. |1999|). Hq and the 



given velocities yield Rvhgo = 21.0 Mpc, which is in good 
agreement with i?viigo = 20.7 Mpc fou nd by Federspiel et al. 
( |1998| ). Note also that Federspiel et al. ( |1998| ) found the same 
value for Hq from relative cluster distances to Virgo. Finally, 
one should recognize that our distance estimate is valid only 
if Virgo is at rest with respect to the FRW-frame. If not, the 
cosmological velocity of Virgo is something more complicated 
than simply the sum of V^q and Vyhto ■ 

4. Results in the direction of Virgo (© < 30°) 

We restrict ourselves to the Einstein-deSitter universe by as- 
signing go = 0.5. As regards Eq. 0we now possess the dis- 
tance d for each galaxy as well as the value k' — 0.606 from 
Eq. O We also have 7 from Eq. O. We need to estimate the 
best value for a. We devised a simple statistical test by finding 
which value or values of a minimize the average | V^bs ^ Vprcd | 
of sample galaxies in the Virgo direction. We found minimum 
values around a = 2.7—3.0. We adopt a — 2.85. The resulting 
systemic velocityH vs. distance diagram is shown in Fig H. The 

' Hq derived is not completely independent because both the SNe 
and the extragalactic PL-relation are calibrated with same Cepheids. 

^ Systemic velocity is a combination of the cosmological velocity 
and the velocity induced by Virgo. This is in our case the observed 
velocity defined in Sect. 4. It could contain also other components, 
but we assume that the Virgo infall dominates. 



IC4182 
NGC3351 
NGC 3368 
NGC 3627 
NGC 4321 
NGC 4414 
NGC 4496A 
NGC 4535 
NGC 4536 
NGC 4548 
NGC 4639 
NGC 4725 



26.4 
26.2 
25.5 
17.3 
4.0 
18.9 
8.4 
4.3 
10.2 
2.4 
3.0 
13.9 



337.0 
640.0 
761.0 
596.0 
1475.0 
692.0 
1574.0 
1821.0 
1641.0 
379.0 
888.0 
1160.0 



0.24 
0.46 
0.52 
0.42 
0.71 
0.80 
0.72 
0.72 
0.65 
0.73 
1.14 
0.58 



345.0 
666.0 
753.0 
698.0 
1611.0 
1098.0 
1482.0 
1642.0 
1261.0 
1710.0 
-409.0 
1025.0 



0.24 
0.47 
0.53 
0.43 
0.72 
0.81 
0.73 
0.73 
0.66 
0.74 
1.15 
0.59 



335.0 
545.0 
728.0 
680.0 
1593.0 
1040.0 
1453.0 
1623.0 
1235.0 
1695.0 
-351.0 
1001.0 



Table 1. Parameters for the 12 galaxies within Q < 30° from 
the adopted centre of the LSC. Column 1: Name, Column 2: 
Angular distance 8, Column 3: Observed velocity V^bs, Col- 
umn 4: Relative distance with i?virgo = 21 Mpc, Column 5: 
Predicted velocity from Model 1, Column 6: Relative distance 
with i?virgo = 20.7 Mpc, and Column 7: Predicted velocity 
from Model 2. 



observed velocities are labelled with circles and the predicted 
values with crosses. The straight line is what one expects from 
Hubble law with our choice of Ho and the curve is the theoret- 
ical TB-solution for 9 = 7° (most of the galaxies have small 
angular distances)^ The relative distances and predicted veloc- 
ities are also tabulated in Table 1 at columns 4 and 5 for this 
Model 1. The galaxies follow the overall TB-pattern, with ve- 
locities steeply increasing when the Virgo center is approached. 
Exceptions are NGC 4639 lying at 1.14 and NGC 4548 at 0.73. 

NGC 4639 has a positive velocity when - according to the 
model - it should be falling into Virgo from the backside. Can 
we explain this strange behaviour? We note it is at a small an- 
gular distance from Virgo (9 = 3.0°). Perhaps NGC 4639 has 
fallen through the centre and still has some of its frontside in- 
fall velocity left. There are, however, no reports of significant 
hydrogen deficiency which one would expect if the galaxy ac- 
tually had travelled through the centre. Maybe this galaxy is a 
genuine member of the Virgo cluster thus having a rather large 
dispersion in velocity. In fact, Federspiel et al. (1998) include 
it to their "fiducial sample" of Virgo cluster galaxies belonging 
to the subgroup "big A" around M87. As a matter of fact our 
chosen centre (l — 284°, b — 74.5°) lies close to the massive 
pair M86/87. 

NGC 4548 is also close to Virgo and has a small angular 
distance (9 = 2.4°). Now, however the situation is reversed. 
This galaxy should be falling into Virgo with a high positive 
velocity but has a small velocity. Paper I classifies this galaxy 
belonging to a region interpreted as a component expanding 
from Virgo (cf. Fig. 8 (Region B) and Table Ic in Paper I). 



^ This curve is simply to guide the eye. The actual comparison is 
made between each observed and predicted point. 
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Fig. 1. The systemic velocity vs. distance for galaxies listed in 
Table 1 for the Model 1 . Circles refer to the observed velocities 
and crosses to the TB-predictions. The straight line is the Hub- 
ble law for Hq = 57 kms~^ Mpc~^ and the curve is the the- 
oretical TB -pattern for 8 = 7°. Note that the most discrepant 
galaxies (NGC 4639 and NGC 4548) are actually closest to the 
Virgo centre in the sky. 



Name 


e 


l^ob. 


-Rmpc 


IC 1613 


163.2 


-62.0 


0.69 


LMC 


107.4 


82.0 


0.05 


NGC 1365 


132.7 


1563.0 


18.20 


NGC 2090 


106.0 


754.0 


11.48 


NGC 224 


126.3 


-13.0 


0.87 


NGC 2541 


63.8 


645.0 


11.59 


NGC 300 


154.1 


125.0 


2.17 


NGC 3031 


61.8 


124.0 


3.37 


NGC 3109 


52.7 


130.0 


1.02 


NGC 3621 


48.4 


436.0 


6.61 


NGC 3198 


43.3 


704.0 


13.68 


NGC 4603 


53.4 


2300.0 


32.51 


NGC 5253 


47.0 


155.0 


3.16 


NGC 5457 


45.7 


361.0 


6.92 


NGC 598 


134.4 


68.0 


0.79 


NGC 6822 


110.7 


8.0 


0.45 


NGC 7331 


125.9 


1115.0 


14.39 


NGC 925 


126.4 


782.0 


8.87 


SEXA 


38.8 


118.0 


1.45 


SEXB 


37.9 


139.0 


1.39 



Table 2. Parameters for the 20 galaxies with & > 30° from 
the adopted centre of the LSC. Column 1: Name, Column 2: 
Angular distance 8, Column 3: Observed velocity K>bs, and 
Column 4; Distance in Mpc. 



Finally, because both galaxies are close to Virgo at small 
angular distances even a slight error in distance determination 
could result in a considerable error in velocities. Hence both 
galaxies might follow the TB-curve but error in distance has 
distorted the figure. 

We also tested the exact values given by Federspiel et 
al. ( |1998| ) with a = 2.85. Now V^f;^^ = 920kms^i and 
Rvhgo = 20.7 Mpc yield with the same infall velocity Hq = 
55 kms^^ Mpc^^. We found k' — 0.641. The behaviour of 
the systemic velocity as a function of distance is quite simi- 
lar to Fig. |l[ This solution is named as Model 2 and given in 
columns 6 and 7 in Table 1 . Due to the paucity of the sample it 
is not possible to decide between these models. 

Though we have only a few points, it is quite remarkable 
how well the dynamical behaviour of galaxies in the direction 
of Virgo is demonstrated by a simple model. That the galax- 
ies follow so well a spherically symmetric model despite the 
observed dumpiness of the Virgo region is promising. Fig. |] 
clearly lends credence to a presumption that the gravitating 
matter is distributed more symmetrically than the luminous 
matter. We have given an affirmative answer to the question 
asked in the introduction. 



5. Results outside Virgo (0 > 30°) 

5.1. On the chosen Hubble constant 

In Paper I the value of the Hubble constant was taken to be 
ifo = 70 kms^^ Mpc^^, consistent with the Tully-Fisher cal- 
ibration adopted at the time and with the used Virgo distances 
(16.5 and 18.4 Mpc giving different infall velocities of LG). In 



the present paper we have fixed Hq from more global consid- 
erations as well as V^q corresponding to i?viigo = 21 Mpc. 
As discussed by ET94 once the velocity of the LG is fixed the 
predicted velocities no longer depend on the value of Hq. It 
is still interesting to see whether the sample galaxies locally 
agree with Hq = 57 kms^^ Mpc^^. This is done in Fig. 0. 
where we have plotted the observed systemic velocities (open 
circles) outside the Virgo region (8 > 30°) as a function of the 
absolute distance in Mpc. The data are given in Table 2. The 
galaxies follow quite well the quiescent flow. In particular, the 
closest-by galaxies (i?gai < 4 Mpc) follow the linear predic- 
tion with surprising accuracy. Note also that galaxies listed in 
Table 1 (galaxies partaking in the Virgo infall) would predict 
Hq = 76 ± 9 km s~^ Mpc^^. This too high a value clearly un- 
derlines the need for correct and adequate kinematic al model 
in Virgo region. 

However, two of the galaxies are clearly discrepant (NGC 
1365 and NGC 4603) and two (NGC 925 and NGC 7331) dis- 
agree. First of all, NGC 4603 is a distant galaxy and the dis- 
tance determination may be biased by the Cepheid incomplete- 



ness effect (Lanoix et al. 1999c). It also does not have a good 



PL-relation in V, which means that one should assign a very 
low weight for it. On the other hand, its velocity differs from 
the quiescent Hubble flow only by 375 kms"^ which actually 
is not very much considering the distance. It belongs to the 
Virgo Southern extension and its velocity could be influenced 
by the Hydra-Centaurus complex. NGC 1365 is a member of 
Fornax, NGC 925 belongs to the NGC 1023 group and NGC 
7331 belongs to a small group near the Local Void. Under these 
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Fig. 2. The observed systemic velocity (open circles) vs. abso- 
lute distance for galaxies outside the Virgo region. The straight 
lines are the predictions for Hq — 57 km s^^ Mpc^^ from SNe 
la and for Hq — 75 kms~^ Mpc^^. Corrected velocities are 
also plotted as crosses. Note that NGC 4603 has a very bad PL 
relation in V and hence one should assign a low weight to it. 
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Fig. 3. The mean heliocentric velocities (filled circles) and ve- 
lociti es corr ected to the centroid of the LG (acc ording to Yahil 
et al. 1977 crosses; according to Richter et al. 1987| : circles) 
for the nearest-by galaxies. 



striking effect. In Fig. g one can see how the galaxies after the 
correction obey the Hubble law all the way down to i?gai = 
0. The correction according to Yahil et al. (1977) are shown 
as crosses and according to Richter et al. (1987) as circles. It 
seems that it is a matter of taste which correction one prefers. 



conditions it is not surprising that these galaxies show devia- 
tions from the Hubble law, except that they all show a tendency 
to have larger velocity than predicted by the Hubble law. 

As can be seen from Fig. |[ these galaxies suggest a much 
shorter distance to Virgo (for them Hq ~ 75 km^^ Mpc^^ 
corresponding to i?viigo = 16 Mpc would be more suitable 
choice for the Hubble constant). This means that without fur- 
ther analysis the data outside the 30° cone do not allow us to 
exclude a shorter distance to Virgo. We examine the effect of a 
short distance to Virgo on the galaxies within the 30° cone in 
an appendix to this paper. 

Finally, we have used so far systemic velocities, i.e. veloci- 
ties not corrected for the virgocentric motions. We tested how a 
simple correction would affect Fig. |[ The corrected velocities Mviriai = 7.5 x lO^''M0i?viigo/16.8 Mpc 
shown as crosses in Fig. g were calculated as 



6. The virial mass of Virgo 

6.1. Prediction from Model 1 



Tully and Shaya (1984) estimated the virial mass of the Virgo 



cluster Following the notation of Paper I the mass within d = 
0.105, which corresponds to = 6° at the distance of Virgo, 
is: 



= 9.38 X IQ^^Mq, 



(17) 



Korr = Kbs+[f (d)H-l'(d)] VI " Sm' Q/<P^Vll 0086,(16) 

where v{d) is solved from Eq. |l^ using parameters for our 
Model 1 and the Hubble velocity at the distance d measured 
from the centre of LSC is t'(d)H = K;osm(l) x d. This correc- 
tion does not resolve the problem. 



5.2. The nearest-by galaxies (Rgai < 4 Mpc) 

We noticed in the previous subsection that nearest-by galaxies 
appear to follow the Hubble law quite well. In this subsection 
we take a closer look at these galaxies. To begin with we note 
that the mean heliocentric velocities Vq (shown as filled circles 
in Fig. 0) have a rather large scatter but on the mean they agree 
with our primary choice of Hq. It interesting that these very 
local galaxies agree with our global value of Hq. Furthermore, 
when Vq's are corrected to the centroid of LG one observes a 



where the latter equality is based on the adopted distance 
Rvhgo — 21 Mpc. With the cosmological parameters (Hq = 
57 kms^^ Mpc~^ and qq = 0.5) and with our Model 1 we find 
using (Eq. 14 of Paper I) 



M{d) = 14.76 X qahl 



Ry 



irgo 



16.8 Mpc 



xd^{l + k'd-°'),{U) 



where Hq = iJo/lOOkms ^ Mpc ^, the mass within d — 
0.105 in terms of A/viriai^ 



Mp.cd = 1.62 X A/v 



rial ■ 



(19) 



The mass deduced agrees with the estimates 1.5-2.0 found by 
Paper I, where it was suspected that the virial mass estimation 
may come from a more concentrated area, which could explain 
the higher value obtained from the TB-solution. 
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Fig. 4. The systemic velocity vs. distance for galaxies listed 
in Table 1 with a — 2.634 (Model 3) predicting exactly 
one Virgo virial mass within a radius d — 0.105 from the 
Virgo. Circles refer to the observed velocities and crosses to 
the TB -predictions. The straight line is the Hubble law for 
Hq = 57 kms^^ Mpc^^ and the curve is the theoretical TB- 
patternforG = 7°. 



6.2. An alternative solution (Model 3) 

We know all the parameters except a needed in Eq. |l8[ This 
led us make a test by looking for such an a which would bring 
about exactly one Virgo virial mass within the radius d = 0. 105 
from the Virgo center We found a = 2.634. The systemic ve- 
locity vs. distance diagram for this model 3 is given in Fig. H. 
Again, NGC 4548 and NGC 4639 show anomalous behaviour 
Comments made earlier are valid. There are three other galax- 
ies showing a relatively large disagreement with the given TB- 
solution. NGC 4535 is close to Virgo at a small angular dis- 
tance. Hence the arguments used for NGC 4639 are accept- 
able also here. Note that in Paper I this galaxy was classified 
as belonging to the region A, where galaxies are falling di- 
rectly into Virgo. Furthermore, Federspiel et al. (1998) classi- 
fies this galaxy as a member of the subgroup B (these galaxies 
lie within 2.4° of M49). This clearly suggests a distortion in 
the velocity. NGC 4414 may belong to the Region B in Paper 
I. This could explain its low velocity. It is also rather close to 
the tangential point so that projection effects can be important. 
For NGC 4536 it is difficult to find any reasonable explanation. 
With these notes taken into account we conclude that model 3 
cannot be excluded. 



ter In the present work we avoided the latter problem up to a 
degree by using photometric distances from the extragalactic 
PL-relation. Such distances are far more accurate than Tully- 
Fisher distances. 

It was quite satisfying to find out that the spherically sym- 
metric TB-model predicts the observed velocity field well. As 



2 a matter of fact, in a recent study Hanski et al. (1999) imple 



mented the TB-model to the mass determination of the Perseus- 
Pisces supercluster. Though the irregular behaviour of the lumi- 
nous matter distribution is even more pronounced in Perseus- 
Pisces than in LSC the mass estimates were reasonable. These 
findings tend to indicate that the gravitating mass is more sym- 
metrically distributed than the luminous matter indicates. 

We found a solution (Model 1) predicting a Virgo clus- 



ter mass within d^ 



■21 



0.105 of AL 



V irgo 



1.62 X M, 



virial 



with il/viriai given by Tully & Shaya (1984). This result is in 
agreement with Paper I, where the difference was suspected 
to mean that the virial mass is estimated from a more con- 
centrated volume. Another plausible explanation is that Virgo 
is flattened. When presuming spherical symmetry under such 
condition more mass is required to induce the same effect on 
the velocities. 

We were also able to find a solution predicting exactly one 
virial mass. Though this model does not agree with observa- 
tions as well as Model 1 the fit is - considering the large uncer- 
tainties involved - acceptable. 

We find these results significant in three ways. The in- 
trinsic behaviour of the sample outside Virgo does not dis- 
agree with the global value of the Hubble constant Hq = 
57 kms^^ Mpc^^. The TB-solutions agree with i?virgo = 
21 Mpc, a value in excellent concordance with i?virgo — 



20.7 Mpc given by Federspiel et al. (1998). It is interesting 
also to be able to predict exactly one Virgo virial mass in the 
Einstein-deSitter universe, because in Paper I the mass predic- 
tions were larger than one and because J7o = 1 is strongly sup- 
ported by theoretical considerations (both Inflation and Grand 
Unification require this value). Discussion on the cosmological 
constant A is postponed to a later phase of our research pro- 
gramme. 

Having said all this one should be cautious. Results pre- 
sented in this paper are based on frontside galaxies only. As 
can be seen from Fig. 4. in Paper I, predicted systemic veloci- 
ties are more sensitive to the model parameters in the backside 
than in the front. In the next paper of this series we intend to 
find a few background galaxies with as reliable distance moduli 
as possible from other photometric distance estimators, in par- 
ticular the Tully-Fisher relation with proper care taken of the 
selection effects. 



7. Conclusions 

Studies of large-scale density enhancements in the Local Uni- 
verse using spiral field galaxies - objects suitable for such a 
project - are discouraged by two factors. The observed distribu- 
tion of spirals is rather irregular and the photometric distances 
based on Tully-Fisher relation are uncertain due to large scat- 
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Appendix A: Is .Rvirgo ~ 16 Mpc? 

In the present paper we fixed the distance to the centre of 
LSC by adopting a iJg from external considerations and by 
fixing the cosmological velocity of the centre of LSC leading 
to i? Virgo = 21 Mpc. In Sect. 5.1 we noted that the galaxies 
outside the 30° cone measured from the centre of LSC do not 
allow us to exclude higher values of Rq. Due to the fixed ve- 
locity these higher values, if really cosmological, necessarily 
lead to a shorter distance to the centre of LSC. 

In this appendix we examine what happens to the TB pat- 
tern within the 30° cone where the TB model is relevant if we 
use i?virgo = 16 Mpc. The result is shown in Fig. |A.1[ The in- 
fall pattern is still clearly present. The agreement between the 
observed points (circles) and the predicted points (crosses) is 
not as good as in Figs, n^ and H. When Fig. A. 1 is carefully com- 
pared with Fig. n^one notes that in Fig. A.l the observed sys- 
temic velocities are systematically smaller than the predicted 
ones. This is because for the shorter Virgo distance galaxies in 
front get closer to Virgo and thus the dynamical influence of 
Virgo should be larger. In Fig. |l] we observe no such system- 
atic effect. It is thus possible to say that our sample is more 
favourable to a long distance scale than to a short scale. 

We also examine the behaviour of the Hubble ratios 
Korr/^gai ^s a function of T4orr (cf- Eq. [I^). In Fig. A.. 2 the 
correction is based on i?viigo = 21 Mpc and in Fig. A.. 3 on 
^Virgo = 16 Mpc. The two values of Hq used in Fig H are 
also given as straight horizontal lines. In both diagrams we ob- 
serve increase in the Hubble ratios as Kon increases starting 
at Vcorr ~ 600 kms~^. In particular, for the shorter Virgo dis- 
tance this effect is quite pronounced. Sandage and Tammann 
have on many occasion stressed that the value of the Hubble 
constant, Hq, should not increase without any clear physical 
reason. In the absence of such the increase in Hq is a warn- 
ing signal that there is a bias present (e.g. the Malmquist bias). 
This is also our tentative interpretation for the increase in Hq. 
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Fig. A.2. The Hubble ratios as a function of the corrected ve- 
locity for i?virgo = 21 Mpc. 
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Fig. A.3. As Fig. A.2 but now i?virgo = 16 Mpc. 



One should also pay attention to the fact that in both cases the 
mean Hubble ratio is ^ 55 kms^^ Mpc^^ for galaxies below 
Fcorr = 600kms-i. 

However, part of the apparent increasing trend in Hq could 
be an artefact of the velocity dispersion. If we overestimate the 
velocities galaxies have a tendency to move upwards and right. 
What is the nature of the suspected bias and its relation to the 



incompleteness bias of Lanoix et al. (1999c) as well as what is 



the influence of the velocity effect are amongst the topics we 
discuss in the next paper in this series. 
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